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We calculate Fermi- surface properties of the Cuprate superconductors within the three- 
band Hubbard model using a cluster expansion for the proper self-energy. The Fermi- 
surface topology is in agreement with angular-resolved photoemission data for dopings 
~ 20%. We discuss possible violations of the Luttinger sum-rule for smaller dopings and 
the role of van-Hove singularities in the density of states of the Zhang-Rice singlets. We 
calculate the shift in the chemical potential upon doping and find quantitative agreement 
with recent experiments. 

PACS. 71.10 Many-body theory - electron states 
PACS. 74.20 Hight-T c - theory 
PACS. 79.60 Photoemission 



1 



Introduction In the last years, the experimental characterization of Fermi-surface 
properties of the high-temperature superconductors has progressed with the help of angle- 
resolved^' 2 '^' 4 '^ and angle-integrat ed( 6 ) photoemission studies. The such determined 
Fermi-surface topology agrees quite well with results from band-structure calculations 
based on the local-density approximation^ 7 ). This measured Fermi-surface topology has 
been taken as an input for theoretical calculations^ 8 ) in the so called "Fermiology" sce- 
nario. Many theorists use, on the other hand, a simplified three-band Hamiltonian (also 
called Emery model^ 9 )), which incorporates strong, local correlations, as exemplified by 
the composite-operator approach^). Determination of the Fermi-surface properties of the 

Emery model is of central interest and has been investigated in a series of recent studies 
(11,12,13,14) 

It is difficult to calculate with numerical methods, like exact diagonalization^ 15 ' 16 ) 
or Quantum Monte-Carlo^ 17 ), the Fermi-surface topology for the three-band model di- 
rectly, due to cluster-size limitations. The numerical methods have been shown, on the 
other side, to be very useful for extracting the energy scales of the local charge- and spin 
excitations^ 15 ). Here we propose that the cluster self-energy provides the missing link be- 
tween such numerical approaches and calculation of momentum-dependent properties, like 
the Fermi-surface topology. 

Method The self-energy of finite clusters with open boundary conditions contains in 
a very precise way information upon the Green's function and self-energy of the extended 
system*^ 18 ). Let us consider the real-space components of the self-energy for both the 

extended system, Ejj(u;), and for a given finite cluster with open boundary conditions, 

(c) ' (c) 

£„• .-(a;). It is straightforward, that all irreducible diagrams contributing to E} '(u>) are 

a subset of the diagrams contributing to Ej j(u;). Therefore the cluster self-energy is a 

systematic approximation to the self-energy of the full system. 

(c) 

Typically one first evaluates the cluster Green's function, G^j(u), in an exact diago- 
nalization study. The cluster self-energy is then obtained by inverting the cluster Dyson's 
equation 

GgM = G<f (uO+£G<f ME<l(uOG<>), (1) 

m,n 

where the i,m,n,j are sites of the cluster and G^\uj) is the non-interacting cluster 
Green's function. 

(c) 

The Green's function of the full Cu02 lattice is determined by inserting the £^ n (o>) 
determined from Eq. (1), into the Dyson's equation of the extended system 

G id {u) = Gl%) + Y,Gl%)J$n(u>)G nd (u,), (2) 

m,n 

where the i,m,n,j run over all sites of the Cu02 lattice and G^-iuj) is the non- 
interacting Green's function of the Cu02 plane. The Green's function, Gij{u), obtained 
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by Eq. (2) can be systematically improved considering larger and larger clusters in Eq. 
(1). The smallest cluster one could possibly consider is the one-site-cluster and Eq. (2) 
together with Eq. (1) would then reduce to the Hubbard-I approximation^). For the 
Emery model we consider the CUO4 cluster for Eq. (1). Within this approximation both 
the important local charge fluctuations and the singlet-triplet splitting are included. 

The passage of the quasi-particle peak through the chemical potential, determined from 
the k-dependent Green's function, G k (uj) = 1/L exp[zk(Rj — Hj)]Gij(ui) determines 
then the Fermi-surface. The density of electrons is given by 

/°° — 1 f d?k 

du( — )ImG(u)n F (u- fx) = 2 ——^Z k , (3) 

with np{oj — ji) and Z k being the Fermi distribution function and the quasiparticle spectral 
weight respectively. The factor of two in Eq. (3) comes from spin-degeneracy and the k- 
sum is over all occupied states. The density of electrons, as given by Eq. (3), is to be 

determined self-consistently( 18 ) as a function of the chemical potential, /i, by demanding 

(c) 

it to be equal to the average density of electrons of the CUO4 cluster, for which E> '(u;) 
is determined via Eq. (1). 

Results For the three-band Hamiltonian we have chosen (in electron notation), fol- 
lowing Eskes and Sawatzky( 15 ), e p — = 5.3eV for the difference in oxygen- and copper 
orbital energies^ 20 ), = 1.3eV and t pp = 0.65eV for the Cu-0 and the 0-0 hopping ma- 
trix element respectively and = 8.8eV for the onsite Coulomb repulsion on the copper 
site. 

From the Green's function of a CUO4 cluster we have extracted the self-energy and used 
it as an approximation for the self-energy of the infinite Cu02 layer, as explained above. 
In Fig.l we illustrate typical results for the angle-integrated spectral function and a doping 
x ~ 20%. For comparison in Fig. 1(c) the density of states of the non-interacting Cu02 
plane (i.e. with = 0) is shown. The bonding-band (mainly copper) is located around 
~ — 9eV, the antibonding band (mainly oxygen) extends from ~ — 3eV to ~ leV and the 
non-bonding oxygen orbitals is in between. For every band the typical 2D logarithmic 
van-Hove singularity is clearly visible, for the non-bonding band, for instance, it is located 
at jj, — 5.3eV = — 3.1eV. In Fig. 1(b) the location and the weight of the peaks of the CUO4- 
cluster Green's function are shown. 

The results for the full Green's function of the extended Cu02 layer are presented in 
Fig. 1(a). Coherently propagating many-body bands have formed out of the cluster energy- 
levels^ 18 ). The chemical potential lies in the band of Zhang-Rice sing separated by 

a charge-transfer gap of about ~ leV from the empty band of copper-d 10 states. The non- 
bonding oxygen orbitals and the band of Zhang-Rice triplets are located -3eV to -5.5eV 
below the Fermi-level. 

Van Hove singularities A van Hove singularity can be observed in the band of 
Zhang- Rice singlets (see Fig. 1(a)). It is tempting to identify this structure with the very 
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sharp feature observed in recent angle-resolved photoemission experiments^ 5 ) on untwinned 
YBa2Cu30g 9 about leV below the Fermi-surface. Let us note that, contrary to predictions 
based on a weak-coupling picture, the van Hove singularity in the density of states does 
not occur at half-filling in our many body calculation. We find that only for relative 
large dopings, namely of about 38%, the Fermi-surface would be located right at the van 
Hove singularity. The absence of a van Hove singularity at the Fermi-edge for dopings 
x < 0.2 is consistent with recent NMR experiments^ 22 ) and in contrast to a band-structure 
calculation based on the local density approximation^ 23 ) . Also we note that these results 
do not support the recently proposed "van Hove scenario" of cuprate supercoductivity^ 24 ). 

Fermi surface In Fig. 2 (a) we present the results for the Fermi-surface as a function 
of doping, x = 0.01, . . .,0.70. Also given are values of the quasiparticle spectral weight, 
Z k , along the (1,1) direction at the Fermi edge. In Fig. 2(b) we present, for comparison, 
the Fermi-surface of the non-interacting CuC"2 layer at various dopings. (0,0) corresponds 
to the T-point. 

We find that the Fermi-surface presented in Fig. 2(a) is in agreement with angular 
resolved photoemission experiments^ 1 ' 3 ) for dopings of about x ~ 20%. The location along 
the (1,0) — (1,1) lines does also agree quantitatively quite well. Along the (0,0) — (1,1) line 
the calculated Fermi edge is located at somewhat too large momenta. We note that the 
quasiparticle spectral weight is reduced by a factor 2 to 5, depending on the doping. For 
small dopings the calculated Fermi surface expands rapidly, due to the reduced spectral 
weight, Z k (compare Eq. (3)) and finally closes completely for x = at the (1,1) point. 
Similar results for the volume of the Fermi function of doping has been found by 

Unger and Fulde^ 12 ) with the projection technique. Note that the data presented in Fig.l 
are for a paramagnetic state and are therefore not in contradiction with the notion^ 13 ' 25 ) 
that in the antiferromagnetic state and dopings of a few percent, the Fermi-surface should 
take the form of "hole-pockets" centered around (n/2,n/2). 

For small interaction strength, in perturbation theory, the Luttinger sum-rule^ 2 ^) is 
valid. It states, that the volume of the Fermi-sea of the system with interaction should 
be equal to that of the non-interacting system. Comparing Fig. 2(b) and Fig. 2 (a) we note 
that the volume of the Fermi sea obtained for = 8.8eV is larger than the one predicted 
by the Luttinger sum-rule (i.e. those obtained for = 0.0). E.g. for x = 20% the 
difference is ~ 14% of the Brillouin zone. The cluster-expansion for the self-energy used 
here is best for large interacting strength and it is therefore not so surprising that the 
results for the Fermi-surface violate the Luttinger sum-rule, which may be expected to 
hold for small interactions, U^. Let us note, that for very large dopings the Fermi-surface 
for both Ud = 8.8eV and = have nearly identical volumes (for x = 70% they 
agree within a percent) , since the doped charge carriers go predominantly into the oxygen 
orbitals, which are influenced only indirectly by Uj. 

Quantitatively the Luttinger sum-rule can be experimentally verified best by compari- 
son with results from band-structure calculations. The experimentally determined Fermi- 
surfaces ( accuracy ( 4 ); 10-20% of the Brillouin zone) agree quite well with band-structure 
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calculations ( 7 ). Since, on the other hand, standard band-structure calculations employ 
only one Slater determinant they necessarily fulfill the Luttinger sum-rule (and yield a 
Z k = l). 

Chemical potential In Fig.3 we present the shift (relative to half-filling^ 27 )) in the 
chemical potential, //, upon doping, x. For comparison the shift in the chemical potential 
for the non-interacting Cu02 plane, [iq, is also shown, they differ by about a factor of 
two. We have included in Fig.3 the data obtained by a recent photoemission study*- 6 ) on 
E^S^Cai^Y^C^Og^. In order to compare with the experimental results we adjusted 
our calculated position of the chemical potential for x = to the experimental value for 
x = 0.01. Note that, eventhough the experimental data for x = 0.01 . . .0.03 correspond to 
an insulator state, while we have a metallic state, excluding these points and adjusting at 
values, e.g., x = 0.08, wouldn't change the curves significantly. 

We also show in Fig.3 the doping dependence of the renormalized Fermi- velocity, 
vp/vpQ ~ m/m*. We find, for about x = 20% doping, an effective mass m*/m of 
about two, in agreement with results from angle-resolved photoemission experiments^ 1 ) 
and quantum Monte-Carlo studies^ 17 ). The effective mass is found to diverge for dopings 
x — > (see Fig.3), as a signal of the metal-insulator transition at half filling. 

Conclusions It is still controversial how the Fermi surface of 2D correlated, trans- 
lational invariant electrons on a lattice changes when going from the metal to the Mott- 
Hubbard insulator regime as a function of electron density. While it seems clear ( 13 ) that 
the first holes doped into the half-filled, antiferromagnetic ordered insulator, occupy states 
with momenta around (±7r/2, ±7r/2), the situation for finite doping concentrations in the 
paramagnetic state is not so clear. Here we proposed a scenario, where the volume of the 
Fermi-sea expands continuously with particle density towards the (tt, it) point at half-filling. 

For moderate dopings, x ~ 0.2%, this scenario is compatible with exact diagonalization 
results of one-band models^ 25 ) and angle-resolved photoemission^ 3 ) experiments on the 
cuprate superconductors. For small dopings, x = 5 — 10% we predict deviations from 
the Luttinger sum-rule, which might be measurable. The calculated shift in chemical 
potential associated with the Fermi-surface topology is in reasonable agreement with recent 
photoemission experiments^ 6 ). 
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Figure captions 

Fig.l (a) A typical result for the angle-integrated spectal weight, for a doping x=0.2. The 
Fermi surface lies in the band of coherently propagating Zhang-Rice singlets, (b) 
The position and the weight of the poles of the CUO4 cluster Green's function used 
for calculating the self-energy, (c) For comparison the density of states of the non- 
interacting (Ufi = 0) CuC>2 plane with the three van Hove singularities for the bonding, 
non-bonding and the anti-bonding band. 

Fig. 2 (a) Results for the position of the Fermi surface in the first Brillouin zone, as a function 
of doping, x, and the quasiparticle renormalization factor, Z, at the Fermi edge. Note 
the strong reduction of Z near the (71- , ir) point, (b) For comparison the Fermi surface 
of the non-interacting (U^ = 0) Cu02 plane, which fulfills the Luttinger sum-rule, is 
shown. 

Fig. 3 Filled circles: Results for the shift in the chemical potential, fx, upon doping, x, relative 
to half-filling. Crosses: Experimental^) results for the shift in chemical potential for 
Bi2Sr2Cai_ x Y x Cu20g +( 5. Open circles: For comparison the shift in (iq of the non- 
interacting, Ud = 0, CuC>2 plane. Filled squares: The renormalized Fermi velocity, 
vp/vpQ ~ m/m* . Note the reduction of about two for x = 0.2. 
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